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Contact effects in devices incorporating strongly-correlated electronic materials are comparatively 
unexplored. We have investigated the electrically-driven phase transition in magnetite (100) thin 
films by four-terminal methods. In the lateral configuration, the channel length is less than 2 /im, 
and voltage-probe wires ~100 nm in width are directly patterned within the channel. Multilead 
measurements quantitatively separate the contributions of each electrode interface and the magnetite 
channel. We demonstrate that on the onset of the transition contact resistances at both source and 
drain electrodes and the resistance of magnetite channel decrease abruptly. Temperature dependent 
electrical measurements below the Verwey temperature indicate thermally activated transport over 
the charge gap. The behavior of the magnetite system at a transition point is consistent with a 
theoretically predicted transition mechanism of charge gap closure by electric field. 

PACS numbers: 71.30.+h,73.50.-h,72.20.Ht 



The complex iron oxide, magnetite, Fe3C>4, is an ex- 
ample of strongly correlated 3d-electron systems^. It 
has been known for decades that bulk magnetite under- 
goes a first-order metal-insulator transition (two-order- 
of-magnitude change in electrical resistivity) at the so- 
called Verwey temperature, Ty ^120 K, accompanied by 
a structural transformation 2 . 

Efforts on magnetite characterization are numerous in 
the seventy years since the discovery of the Verwey tran- 
sition, including thorough investigations of its electri- 
cal properties^— supported by theoretical calculations of 
electronic structure^—. Recent advances in nanofabrica- 
tion and film growth allow electrical characterization at 
previously inaccessible scales, leading to the recent dis- 
covery of an electric field driven transition: Magnetite 
films or nanoparticles below Ty experience a transition 
from an insulating state to a state with much lower resis- 
tance upon application of a sufficiently high voltage^— . 
The switching voltage scales linearly with the channel 
length suggesting an electric-field driven transition. 

The key point of these experiments was an examina- 
tion of magnetite films or nanoparticles between two elec- 
trodes separated by only several hundreds of nanometers 
or less. In this configuration the electric field needed to 
drive the transition was accessible at relatively low volt- 
ages, thus preventing both excessive heating and dam- 
aging of the sample. We proved the observed switching 
not to be an artifact of heatin g 11 : 13 , in contrast to previ- 
ously observed transitions in magnetite driven by Joule 
heating of the samples above Ty under bia o 14 : 15 . 

The downside of such small channel length experi- 
ments is an unavoidable, dominant contribution of the 
contacts, which prevents direct insight into the proper- 
ties of magnetite before and after transition. By fitting 
our data for two-terminal devices with different chan- 



nel lengths it was demonstrated that contact resistance 
of Au/magnctitc interfaces comprises more than 70 % of 
the total resistance^. Upon testing several different con- 
tact metals (Au, Pt, Cu, Fe and Al), copper showed the 
lowest contact resistance with magnetite filmH. Even 
with a Cu contacting layer, however, the contribution of 
the contacts to the total two-terminal device resistance 
cannot be neglected. 

One of the most effective ways to differentiate between 
bulk and interface effect is to make multilead measure- 
ments. To date no such experiments have been performed 
to study the recently discovered electrically driven transi- 
tion in magnetite. In this paper we perform four-terminal 
experiments in a lateral electrode configuration using 
magnetite thin films. The channel length is less than 
2 /nm and voltage-probe wires ~100 nm in width are 
directly inserted into the channel. These multilead ex- 
periments quantitatively and unambigously separate the 
role of each interface and the magnetite channel. For 
the first time we study the changes in contact and chan- 
nel resistance contributions at the onset of electric-field- 
driven transition in magnetite. Results indicate that at 
the transition point both contact resistances and the resis- 
tance of the magnetite channel decrease abruptly. By do- 
ing temperature-dependent electrical measurements be- 
low Ty we trace the thermally activated transport over 
the charge gap in magnetite and provide an insight into 
the transition mechanism in this system. 

The Fe 3 4 (100) thin films (thickness: 50-100nm) used 
in the present study were grown on (100) oriented MgO 
single crystal substrates using oxygen plasma assisted 
molecular beam epitaxy system (DCA MBE M600) with 
a base pressure 2xl0 -10 Torr. The substrates were 
cleaned in-situ at 873 K in 5xl0~ 6 Torr oxygen for two 
hours. 
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FIG. 1. (color online) (a) SEM image of the device for four- 
probe measurements showing source and drain leads and two 
pairs of voltage probes within the channel, (b) Colored SEM 
image demonstrating electrical contacts to /im-size Au pads 
with further In soldering to attach Au wires, (c) Schematics 
of electrical circuit of four-probe measurements. Letters S and 
D denote source and drain contact, respectively. Contacts are 
made of 6 nm Cu adhesion layer (reddish) and 10-20 nm cover 
layer of Au (yellow), (d) Temperature dependence of the low- 
bias resistance of magnetite channel (-Rdev) and correspond- 
ing contact resistances (-Rc) at source and drain electrodes. 

Reflection high energy electron diffraction, RHEED, 
(STAIB Instruments) was used to monitor the growth 
mode and growth rate (0.3 A/s). Room temperature 
Raman spectroscopy (performed in the backscattcring 
configuration using Rainshaw 1000 Micro Raman sys- 
tem), High resolution X-Ray diffraction measurements 
using a multi-crystal high-resolution X-ray diffractome- 
ter (HRXRD, Bede-Dl, Bede, UK), and low temperature 
four probe resistance measurements were performed to 
establish the crystal structure and stoichiometry of the 
Fe 3 4 phased. 

Devices for two- and four-terminal measurements were 
prepared by electron beam lithography (Jeol-6500 SEM). 
A channel length of 400 nm - 1.9 /im is defined by two 
10 /xm wide source and drain leads. One or two pairs 
of voltage-probe leads were directly patterned within the 
channel (Fig. QJi) . The contacts were fabricated by the 
electron-beam deposition. In a typical experiment 6 nm 
of Cu (the best adhesion to magnetite film and lowest 
contact resitancc out of studied contact metals^) and 
10-20 nm cover layer of Au were deposited. The leads 
were connected to micrometer-size pads (300 x 300 /Ltm) 
to which Au wires are attached by In soldering (Fig. [TJa) 
and then connected to external contacts of the puck. The 
puck was placed into the chamber of a Quantum Design 
Physical Property Measurement System (PPMS model 
6000) for variable temperature measurements (300 K - 
80 K). The lower temperature bound of these measure- 
ments is limited by the increasing switching voltages and 
concerns about device damage as T is further decreased. 

Electrical characterization of the devices was per- 
formed by standard four-terminal methods using a 



semiconductor parameter analyzer (HP 4155A). The 
schematic of device electrical connections is presented in 
Fig. [Tb. The voltage, lout! is applied to the source lead 
with the drain grounded, and current flowing through 
the channel is recorded. The pair of voltage probes, di- 
rectly inserted into the channel between source and drain 
leads, senses voltages V\ and V^. A voltage drop in the 
channel without the contact contribution is then calcu- 
lated as AV = Vi — Va- Only one pair of voltage probes 
(either left or right, see Fig.QJi) is active in a given mea- 
surement, with the second pair being intact. Having two 
pairs of voltage probes allows two independent sets of 
measurements (one for each pair of voltage probes) in a 
given channel, to verify data consistency in these four- 
terminal devices. As was demonstrated in our previous 
paper, sweeping voltage in a continuous staircase mode 
leads to the overheating of the sample and appearance 
of a hysteresis in forward and reverse bias sweeps^. To 
minimize Joule heating of the channel, the voltage was 
always swept in a pulsed regime with the shortest avail- 
able pulse duration (500 [is) and with pulse period > 5 
ms. This pulse measurement procedure greatly reduced 
apparent hysteresis in the transition as a function of bias 
sweepi^. 

At any V ou t, voltage first drops at source 
electrode/Fe304 interface (Vc (source)). We assume 
that at low source-drain biases the contact interface 
contributions are dominated by an Ohmic contribution, 
Rc (source) = Vc (source)//. Then, in the assumption of 
a homogeneous film (medium) between the electrodes, 
voltage linearly drops across the channel, and two values 
are recorded at the two locations of the voltage probes. 
The remaining potential drop to zero volts (grounded 
drain electrode) occurs at Fe304/drain electrode inter- 
face (Vc(drain)— ► i?c (drain)). Conventionally, the total 
device may be represented as a voltage drop over three 
resistors in series, Rc (source), -Rdev and Rc (drain). By 
knowing the geometrical characteristics of our devices 
from SEM images (i.e., t\, A£, and £ 2 , see Fig. [TJ2) 
we can calculate the values of all three voltage drops 
and, by dividing over measured current, corresponding 
resistances. 

An example of the temperature dependence of -RdeVj 
Rc (source) and Rc (drain), calculated this way at V out 
— lOOmV, is presented in Figure [TJi in the temperature 
range around Ty The Verwey temperature is inferred 
for each device as an inflection point in i?DEv(T) depen- 
dence; and for the various devices studied here Ty range 
from 100 K to 110 K. For source and drain electrodes 
made of the same metal (Cu in this case), Rc (source) w 
Rc (drain) , besides, -Rdev and Rc have nearly identical 
temperature dependence (Fig. [T}1 ). 

Let us consider the relative contributions of contacts 
and magnetite channel to the total voltage drop. Our as- 
sumption of linearity in the channel conduction presumes 
that Rdev linearly scales with the channel length, L, 
while Rc (source) and Rc (drain) should remain indepen- 
dent of L. This is supported experimentally. To demon- 
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FIG. 2. (color) (a) Examples of two-terminal resistance de- 
pendence on the channel length at three different tempera- 
tures (85 K, 90 K and 95 K) with corresponding linear fits. 
Extrapolation to zero channel length gives the total contact 
resistance, Rc (total), at each temperature, (b) and (c) show 
SEM images of devices with different channel lengths, (d) 
and (e) plot calculated values of device resistances, Rdev, 
and contact resistances, Rc (source) and Rc (drain), around 
zero source voltage for devices in (b) and (c), respectively. 

strate this, we made a set of devices on the same piece 
of magnetite film with different channel lengths while 
all other geometrical parameters (film thickness and the 
width of source and drain leads) remained exactly the 
same for all devices. Two rcprcsentivc SEM images of 
such devices with L=l.l fixn and L=1.9 /xm are shown 
in Fig. [2b and c, respectively. The total (two-terminal) 
resistance at each temperature linearly depends on the 
channel length as demonstrated in Fig. [3^ at several tem- 
peratures (85 K, 90 K and 95 K). 

The calculations of contact resistances based on l\ , A£ 
and £2 for each device show that at each temperature 
Rc (source) and Rc (drain) are equal to each other and 
are the same for devices with different lengths, L. It is 
worth mentioning that Rc (source) and Rc (drain) do not 
change upon switching the grounds, i.e., exchanging the 
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FIG. 3. (color online) Examples of I-V and corresponding 
I-AV curves at 80 K and 85 K at the voltage ranges above 
switching voltage at each temperature. Arrows indicate the 
direction of the voltage sweep. The pulse sweep parameters 
are those to demonstrate the small hysteresis in forward and 
reverse voltage sweeps. Top inset shows I-AV curves at the 
same temperatures, but at the voltage range below switching 
voltage. Bottom inset zooms in to I-AV curve at 80 K around 
transition point to demonstrate the discontinuity in measured 
AV" value. 

place of injecting and grounded electrodes. Calculated 
resistances in Fig. [2}i and e represent the data for the 
devices shown in Fig. [2b and [2fc, respectively. While 
Rc (source) and Rc (drain) remain independent of chan- 
nel length, i?DEV increases as L increases which is ob- 
vious from comparison of Fig. [2ji and Fig. [2^. All three 
resistances increase significantly with decreasing the tem- 
perature (compare data at 85 K, 90 K and 95 K in Fig.[^H, 
e) , as will be discussed below in detail. The contact resis- 
tance, Rc (source) + Rc (drain), contributes from 20% to 
13% of the total two-terminal R for devices with channel 
lengths ranging from 1 /im to 2 /xm. 

The increase in total two-terminal resistance with 
channel length (Fig. [2^) is caused by the increased con- 
tribution of i?DEV i n longer devices. Moreover, the ex- 
trapolation of a two-terminal R vs L linear fit to zero 
channel length gives a resistance value very close to the 
sum of calculated Rc (source) and Rc (drain) at each tem- 
perature. The latter proves the consistency of our calcu- 
lations and independence of contact resistances on the 
channel length within L range investigated in this work. 

At temperatures below Ty ~ 105-ftT, the current- 
voltage characteristics, I-V, show Ohmic behavior at low 
source voltage range (< IV), while start to exhibit non- 
linearities at higher voltages, symmetrical for positive 
and negative source voltages. Examples at two selected 
temperatures (80 K and 85 K) are shown in the top inset 
of Figure [3l Upon further increasing source voltage I-V 
curves show a sharp jump in current (Fig. [3|) as soon as 
the source voltage reaches a critical switching value, Vsw 1 
at a certain temperature as described in detail i n 11 ! 13 . 
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FIG. 4. (color online) (a) Schematic diagram of voltage dis- 
tribution along the channel at a transition point (blue open 
squares) and right after transition (red closed squares). Vc 
denotes the voltage drop at the interfaces, (b) A fragment of 
I-V curve at 85 K in the vicinity of the transition. Blue open 
square marks the transition point and red closed square shows 
a point right after transition; blue and red squares in (a) cor- 
resond to the voltage distribution over the channel at these 
points, (c) and (d) are the voltage dependences of Rc (source), 
Rc (drain) and Rdev, respectively, demonstrating the abrupt 
decreases (jumps) in all three resistances at the transition 
point. 

This is a transition from high resistance (Off) state to a 
state with much lower resistance (On) state. Note the 
two-order of magnitude difference in current after transi- 
tion by comparing I-V curves before (top inset) and after 
transition in Fig. [3] 

Corresponding / vs Ay = Vi—Vi plots (Fig. [3]) have 
much lower switching AVsw values and reveal at a tran- 
sition point not only a discontinuity in current, but also 
in the measured AV value, which decreases in absolute 
magnitude (Fig. [3] bottom inset). Since Ay = V% — V2, 
in general, reflects properties of magnetite channel with- 
out contact effects, the discontinuity (jump) in Ay at a 
transition point can be explained as a sudden decrease in 
device resistance, -Rdev- 

Now let us turn to the quantitative description of 
contact effects at the onset of the field-induced transi- 
tion. Calculations of Rc (source), Rdev and Rc (drain) 
show that at a transition point the voltage drops at the 
contacts, Vc (source) and Vc (drain), increase in abso- 
lute value, while AV decreases. Compare the blue open 
squares (at the transition point) and red closed squares 
(at the next point after transition) in Fig. [4ji, which de- 
pict voltage distribution over the channel length. From 
this sketch the decrease in AV value is also clearly visible. 

Although Vc (source) and Vc (drain) increase, due to 
the overall increased current, R c (source) and Rc (drain) 
actually decrease upon passing through the transition 
point. Fig.[4}),c, and d explicitly demonstrates these de- 
creases (jumps) in -Rdev and Rc at a transition point, 
denoted further as Rj^f and R c ump . For source and 
drain electrodes made of Cu, Rc (source) and Rc (drain) 
jumps at the transition point are equal to each other and 
remain unchanged in the experiments on exchanging the 
grounds. 

Note that at the transition point both device and 
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FIG. 5. (color online) Temperature dependences of (a) the 
jumps in Rc (source), Rc (drain) and in -Rdev at a tran- 
sition point; each point represents an average over 8 in- 
dependent measurements (left/right voltage pairs, different 
grounds, positive/negative switching voltages), standard de- 
viation is within the symbol size (b) jumps in AV (squares) 
at a transition point and its exponential fit (solid line) (c) 
Vsw and AVsw with corresponding exponential fits and (d) 
conductance (1/Rdev) of magnetite channel in Arrhenius co- 
ordinates and its linear fit. 

source and drain contact resistances decrease abruptly. 
This behavior is distinct from the one for other systems 
exhibiting voltage-driven transitions (such as mangan- 
ites and doped SrTiOa). For these systems the lead- 
ing role of oxygen vacancies drift under applied bias 
was demonstrate d 18 ' 19 , and source and drain contact re- 



sistances show variations of opposite sigr>2£i; 



while 



source contact resistance increases, the drain contact re- 
sistance decreases. Since this is inconsistent with our 
observations, magnetite clearly must exhibit a different 
switching mechanism. 

As temperature decreases, the jumps in contact resis- 
tance, R c ump , and device resistance, R^^y, remain nega- 
tive, but increase in absolute magnitude (Fig.[SJi). Upon 
approaching the temperature when switching is not ob- 
served (T ~ T\ Ul13 ), R 3 c umv and R J jjev approach zero. 
The jump in AV is also temperature dependent and its 
magnitude exponentially decays with the temperature, 
approaching zero at Ty. Vsw and AVsw depend on 
temperature exponentially as well, as demonstrated in 
Fig. [Sk. 

To explain the temperature dependence of the parame- 
ters related to the observed transition, i.e., jumps in cur- 
rent, AV, R c and -Rdev (Fig- El- we should review the 
properties of magnetite below Verwey temperature, since 
the transition is only observed below Ty. We will con- 
sider the Verwey transition physics in magnetite from the 
electronic structure point of view. While magnetite has 
strong electron correlations, it is believed that a band- 
type description of its electronic structure is a reason- 
able approximation^, with transport being dominated by 
low-lying electronic states near an effective Fermi energy, 
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Ep • 

The electronic structure of magnetite has been probed 
extensively by photoclcctron and scanning tunneling 
spectroscopies and band-structure was calculated using 
different method s 10 ' 21 ' 22 . Above Ty there is a finite (non- 
zero) density of states (DOS) around Ep, which is dom- 
inated by Fe 3d states of the -B-site sublattice of cubic 
structure. Below Ty the DOS near Ep exhibits a clear 
gap, causing two orders of magnitude increase in resis- 
tivity at Ty2L22.. 

Recently, it has been theoretically predicted that a gap 
to charged excitations (charge gap) in correlated insu- 
lators can be closed by applying external electric field, 
resulting in field-induced metal-insulator transitio n 24 ' 25 . 
Several systems exhibit this behavior, for example, 
charge-ordered state of complex manganese oxide s 26 ' 27 , 
as well as ID cuprates^ and 2D nickelates 2 ^. 

Our data on magnetite seem to be another experimen- 
tal observation consistent with this sort of gap closure by 
electric field. The absence of hysteresis in forward and 
reverse bias sweeps, meaning that metallic state persists 
only if applied voltage (electric field) exceeds a critical 
value, \V\ > \Vsw\ , is expected from this mechanism 
and is indeed observed in Fe3 04 (see Reb^ 3 - for details). 
The jump in |-Rdev| at the transition point is a natural 
consequence of the gap closure. The increase in the abso- 
lute value of -Rqev as ^ decreases (Fig. [5h) is also easily 
explained, as at lower temperatures there is a transition 
from more insulating state to the same metallic state 
with zero-size gap. Accompanying drops of contact re- 
sistances, R^ mv ', are direct consequences of gap closure 
and, thus, the change in the position of magnetite Ep 
relative to the E-p of the contact metal. 

The exponential dependences of the above parameters, 
particularly that shown in Fig. [5]i, imply thermally ac- 
tivated transport below Ty. Indeed, plotting the inverse 



of .Rdev ( m the Ohmic regime near zero bias) in Ar- 
rhenius coordinates gives a straight line in a given T 
range (FigJSJi). The activation energy, E a , inferred from 
these data on the magnetite channel (device) lies in 85- 
89 meV range for several devices. These values match 
well with the size of the gap below Ty, inferred from 
photoemissio n 21 ' 23 and optica l 30 ' 31 spectroscopies data. 
This suggests that transport below Ty involves charge 
carriers thermally activated over the gap. 

In conclusion, by doing four-terminal experiments at 
magnetite thin films below Ty we quantitatively sepa- 
rate the contributions of each electrode and the mag- 
netite channel before and after the electric field driven 
transition. For devices of increasing channel lengths we 
demonstrate the increase in total resistance to be caused 
by increased contribution of the magnetite channel, while 
contact resistances are unchanged for all channel lengths 
within 1 to 2 fj,m range. At all temperatures the tran- 
sition is observed (T < Ty), contact resistances of both 
source and drain electrodes and the resistance of mag- 
netite channel decrease abruptly at the transition point. 
This behavior is consistent with the mechanism of charge 
gap closure by electric field predicted in theor y 24 ' 25 . To 
further explore the field-driven switching mechanism in 
magnetite, the effect of contact metals with different 
work functions is currently under study. In the frame- 
work of the charge gap closure mechanism, the magnitude 
of contact resistance jumps at a transition point, i?p" mp , 
are expected to be dependent on the work function of 
the contact metal according to the relative alignment of 
metal Fermi level and effective Fermi level of magnetite. 
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